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Abstract 

We recently considered decays in the framework of chiral per- 
turbation theory based on the effective Lagrangian including mesons, 
photons, and leptons. There, we published analytic one-loop-level ex- 
pressions for form factors / and g corresponding to the mixed process, 
K° — > 7T°ir~£ + i'£. We propose here a possible splitting between strong 
and electromagnetic parts allowing analytic (and numerical) evalua- 
tion of Isospin breaking corrections. The latter are sensitive to the 
infrared divergence subtraction scheme and are sizeable near the tttt 
production threshold. Our results should be used for the extraction 
of the P-wave iso- vector tttt phase shift from the outgoing data of the 
currently running KTeV experiment at FNAL. 
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1 Introduction 

Every time that a kaon decays into a couple of pions and a lepton-neutrino 
pair, a tttt scattering occurs in the final state. Whenever a pion scatters on its 
twin, it offers to us an additional opportunity to scrutinize the fundamental 
state of strong interaction (see [I] for references). Let 5f be the phase of a 
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two-pion state of angular momentum / and Isospin / and consider the 
decay process, 

K{p) — > n( Pl )n(p 2 )l + (p e )i; e (p u ) , (1) 

where the lepton, £, is either a muon, /i, or an electron, e, and v stands for 
the corresponding neutrino. In the Isospin limit, the decay amplitude, A, for 
process can be parameterized in terms of three vectorial (F, G, and R) 
and one anomalous (H) form factors, 

G 

A = i -j= KfiipOlvtt- - l 5 ) v (Pe) x 

[(pi + P2YF + (p! - pz fG + [p t + p v yR] 
e^(n + pMPi+P2) P (pi-P2),H} , (2) 



M K ± 
1 



where V us denotes the Cabibbo-Kobayashi-Maskawa flavor-mixing matrix el- 
ement and Gp is the so-called Fermi coupling constant. Note that form 
factors are made dimensionless by inserting the normalizations, M^} ± and 
M^|. The fact that we have used the charged kaon mass is a purely con- 
ventional matter and corresponds to the choice of defining the Isospin limit 
in terms of charged masses. In the following, we will be interested only in 
two form factors, F and G, and denote by (F, G) + ~ and (F, G)°~ the ones 
corresponding to the physical processes, 

K+(p) — n + ( Pl )7r-(p 2 )£ + (p e )u e (p u ), (3) 

and 

K°(p) — 7r°(p 1 )7r-(p 2 )£+(p £ )^(p I/ ), (4) 

respectively. 

Form factors are analytic functions of three independent Lorentz invari- 
ants, 

Stt = (P1+P2) 2 , S e = {pt+PuY, (5) 

and the angle 9^ formed by p 1 , in the dipion rest frame, and the line of flight 
of the dipion as defined in the kaon rest frame PIE]- It has been shown in j3] 
that, in the experimentally relevant region, the partial wave expansion, 

F + - = (fsM + f t s e )e i5 °°^ + fpXY cos 9„e iS ^\ (6) 
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G + - = (gp + g'ps^ + ges e )e is ^ + g D XY cos 9^ iS °^ , (7) 
is proving sufficient to parameterize form factors. In the preceding, 

X = ^ A 1/2 (s„-, si, M 2 K± ) , y = -\^ 2 (s n ,M 2 ± ,M 2 ± ), (8) 

with, 

X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz , (9) 

the usual Kallen function. Note the linear dependence of the first term in 
the partial wave expansion of form factors on Sf. Isospin symmetry, Bose 
symmetry, and the AI = 1/2 rule lead to, 

F°~ = V2f P XY cos e n e iS ^\ (10) 
G°- = V2{gp + g pS „ + g iSi )e i51 ^. (11) 

It follows that Km decay of the neutral kaon is dominated by P waves. 
Therefore, a precise measurement of form factors for the decay in question 
would allow an accurate determination of the P-wave iso-vector nn phase 
shift. 

The currently running KTeV experiment [S] aims at measuring form fac- 
tors for decay of the neutral kaon with an accuracy 3 times better than 
the one offered by previous measurement jHllZj. The outgoing data on form 
factors contain, besides strong interaction contribution, a contribution com- 
ing from electroweak interaction. The latter breaks Isospin symmetry and is 
expected to be sizeable near nir production threshold & . In order to extract 
mi scattering parameters from the KTeV measurement, Isospin breaking cor- 
rection to form factors should therefore be under control. In this direction, 
we recently published analytic expressions for F°~ and G°~ form factors cal- 
culated at one-loop level in the framework of chiral perturbation theory based 
on the effective Lagrangian including mesons, photons, and leptons [Tj. In 
the present work, we will split analytically Isospin limit and Isospin breaking 
part in form factors allowing a first evaluation of Isospin breaking effects in 
Ku decays. 

2 Kinematical variables 

In the following, we shall consider process (jlj and use, unless mentioned, 
notations of reference p. In the presence of Isospin breaking, the decay 
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amplitude for process (j3J) can be written as follows by Lorentz covariance, 



A 




\- [(pi + P2Yf~ + (pi - pa) V~ + ( P i + PuYr Q ~] 7m 



i 



e^{p e + Pv ) u ( Pl + p 2 ) p (p 1 - p 2 ) a ti 



o 



1 



[l^lv] PlP2 T \ V (Pl) 



2Ml 



The quantities /, g, r, and h, will be called the corrected form factors 
since their Isospin limits are nothing else than the form factors, F, G, R, 
and H, respectively. The tensorial form factor T is purely Isospin breaking 
and does not contribute to the mixed process at leading chiral order. The 
corrected form factors as well as the tensorial one are analytic functions of 
five independent Lorentz invariants, s n , se, 6 n , 6?, and <fi. 8e is the angle 
formed by p e , in the dilepton rest frame, and the line of flight of the dilepton 
as defined in the kaon rest frame. <fi is the angle between the normals to 
the planes defined in the kaon rest frame by the pion pair and the lepton 
pair, respectively. Let us denote by 5F and 5G the next-to-leading order 
corrections to the F°~ and G°~ form factors, respectively, 



The analytic expressions for 5F and 5G were given in pQ. We shall dis- 
tinguish between photonic and non photonic contributions to SF and SG. 
The photonic contribution comes from those Feynman diagrams with a vir- 
tual photon exchanged between two meson legs or one meson leg and a pure 
strong vertex. Obviously, this contribution is proportional to e 2 , where e is 
the electric charge, and depends in general on the five independent kinemat- 
ical variables, s n , si, 9^, 0£, and through Lorentz invariants like (p 2 +Pi) 2 , 
say. The non photonic contribution comes from diagrams having similar 
topology as the ones in the pure strong theory with Isospin breaking allowed 
in propagators and vertices. This contribution generates Isospin breaking 
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terms proportional to the rate of SU{2) to SU(3) breaking, 



. V3m d -m u A . 1 

e = — , m = - [m u + m d ) , (12) 

4 m s — m 2 

and to mass square difference between charged and neutral mesons, 

A n = M 2 ± - M 2 = 2Z e 2 F 2 + 0(p 4 ) , (13) 
A K = M 2 K± - M 2 K0 = 2Z e 2 F 2 - B (m d - m u ) + 0(p 4 ) , (14) 

or equivalently, {m d — m u )/(m s — rh), Z e 2 , and m d — m u . The kinematical 
dependence is on three Lorentz invariants, (pi +P2) 2 , (p — Pi) 2 , and (p — p 2 ) 2 
which represent respectively the dipion mass square, the exchange energy 
between the kaon and the neutral pion, and that between the kaon and the 
charged pion. In terms of independent kinematical variables, the preceding 
scalars are functions of s n , si, and cos^. 



2.1 The photonic contribution 

A generic term in the photonic contribution can be, 

photonic contribution = e 2 ((P2 + Pe) 2 , • • •) , (15) 

i 

where £j is an arbitrary loop integral function of (p 2 + Pe) 2 - To the order we 
are working, that is, to leading order in Isospin breaking, the power counting 
scheme we use dictates the following on-shell conditions to be used in the 
argument of £j, 

p 2 = M 2 K = B {m s + rh) , p\ = p\ = M 2 = 2B m . (16) 
Therefore, (p2+Pe) 2 in (fTB]) should be replaced by the following expression [I], 
Ml + m 2 




If m 2 \ ( 4M 2 ^ 1/2 



4 V Si J V s 



[M\ — si — s^) cos 6 n cos 6e 



m 2 \ / 4M 2 ^ 1/2 



+ - (l - J (l - — (s^) 1/2 sin^sin# £ cos0. 

From the foregoing, it is clear that for sg = mf the photonic contribution does 
not depend neither on 6t nor on 0. In order to reduce the complexity of the 
study and allow the treatment of photonic and non photonic contributions 
to A°~ on an equal footing, we will assume that, 

st = mj, (17) 

and use for (p 2 + Pe) 2 in (|T5Jl the following expression, 

(P2 + Pe) 2 = l -(M 2 K + 2Ml + m 2 f -s^) 
1 / 4M 2 \ 1/2 , 

jfl-^J A 1 / 2 ^, m 2 ,M|) cos ^. (18) 

It follows that (|15[1 can be written as, 

photonic contribution = e 2 <?(s w ) + e 2 ^(s^) cos9 n , (19) 
where q and $ are analytic functions of s n . 



2.2 The non photonic contribution 

In order to split strong and electromagnetic terms in the non photonic con- 
tribution, one has to expand the exchange energies, (p — pi) 2 and (p — P2) 2 , 
in powers of the fine structure constant, a, and m d — m u . To this end, we 
shall first express these scalars in terms of and cosO^ for se = m 2 and in 
the presence of Isospin breaking. From pQ , 

(p- Pl ) 2 = M 2 KO + M 2 o 

~ ^[(M 2 K0 -m 2 + S7r )( S7V + M 2 -M 2 ± ) 

+ A 1 / 2 ( S7r ,m 2 ,M| )A 1 / 2 (^,M 2 o,M 2 ± )cos^] , (20) 
(P-P2) 2 = M 2 K0 +M 2 ± 
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— [(M 2 - m\ + - M 2 , + Ml 



- A 1/2 (s 7r ,m, 2 ,M| : o)A 1/2 (^,M2 ,M2 ± )cos^] . (21) 

Let us denote by and u n the Isospin limits of the preceding Lorentz scalars, 
1 



(M 2 K± + 2M 2 ± + m 2 - s 7 



in 



\ 1/2 (s n ,m 2 e ,M^ ± )cos6 w 



ii, 



- (M 2 K± + 2M"U + m 2 - sj 



4M 2 ± 



1/2 



A 1/2 (s 7r , m 2 t ,M 2 K± )cos9^ 



For completeness, it is convenient to note the following proposition, 



cos^ = => U 



it, 



- (M 2 K± + 2M 2 ± + m 2 - s 7 



Using the replacements, 



M 2 



M 2 ± 



M 



K o — > M K ± — A K , 



and expanding (j2"Uj) and (}2*T]) to first order in A n and A K , we obtain, 



(22) 



(23) 



(24) 



(25) 



{P-Pif = \{M 2 K± + 2Ml ± +m 2 i -s TT ) 



+ 
+ 



^(M|-m 2 -^)A^-iA A - 



4-- 

1 



2s w 
1 



1 - 



4M 2 ± 
4M 2 



1/2 



A 1/2 (^,m 2 ,M^ ± ; 



-1/2 



M 1 " „. 



$7V 



A 1/2 (^,m 2 ,M|)A 7 



(M 2 - m 2 - Sw )A- 1 / 2 (^, m 2 , M| r )A K 



cos^, (26) 
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]-(M 2 K± + 2M 2 ± +m 2 -s n ) 
^(M 2 K -m 2 + s^)A 7r --A K 



+ 





[M 2 K -m\- s 1 ,)\' xl2 {s^vn\,M 2 K )A K cos 9. 



(27) 



Note that terms of order O^A^Ax) are forbidden by our power counting 
scheme since they are first order in Isospin breaking. Although equations 
(I26|) and ()27|) are simple to derive, their utility is of great importance to the 
present study. In fact, the involved expansion could be generalized to any 
Ki4 observable as we will see below. 

2.3 Splitting strong and electromagnetic interactions 

The first step in our program consists on injecting equations (I26|) and (J27|) 
in the non photonic contribution to the decay amplitude A°~ . Then, we 
expand once more to first order in A^ and Ak dropping out terms of order 
0(A 7T Ak)- As a result, form factors for Ku decay of the neutral kaon can 
be written in the following compact form which shows explicitly the splitting 
between strong and electromagnetic interactions, 



x 



,o- 



(stt, - Pi) 2 , (p - P2) 2 , (P2 + Pi) 2 , ■ ■ •) 
^ [S xg + U x (s n ) + V x (s T ) cos^] , 



x = f, 9, 



(28) 







where 



X 



W s x + W^Ak + W X K A K 

+W x ie 2 + W e x ^=, W = U , V , 



(29) 
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are analytic functions of s ff . If one makes the following substitutions, 

A u — > 2Z e 2 F 2 , (30) 
A K — 2Z e 2 F 2 - i| (M 2 - M 2 ) , (31) 

then, equations (J2*8|l and (}2T?j) read, 



€ 



w* = w: + wy + w*-= } (32) 



m d -m u /g 

= U^+2Z F 2 (Vy;+l4^), (33) 

m,j-m„ 



W^-m. = ^-4(M|-M 2 )^. (34) 



The aim of the present work is to determine the U functions corresponding 
to / and g form factors for Ku decay of the neutral kaon. 



3 The photonic contribution 



From now on, we will work under proposition ()24|) keeping in mind that, in 
the Isospin breaking contribution, the power counting dictates the following, 

Isospin breaking — > t n = u T = - (Ml + 2M 2 + m 2 — s n ) . (35) 

Taking the photonic contribution from P, applying assumption ()17j) . and 
performing the preceding expansion, it is easy at a first sight to derive U e 2. 
The problem is that, in practice, one encounters loop integrals with vanishing 
Gramm determinant when reducing vector and tensor integrals to scalar 
ones p. After a long and tedious calculation one obtains, 

U[ 2 = l(-6K 3 + 3K 4 + 2K 5 + 2K 6 -QX 1 ) 
3 

2 M 2 

[6K 3 - 3K 4 - 2K 5 - 2K 6 + 2K 9 + 2K W ) 



3 Ml- M* 
+ £?(M 2 ,0,M 2 ) 

- 1 [{M 2 K - m 2 - srf - Am]Ml] X~\u, m 2 , M 2 ) 
+ - A mjMl(s n - AMl){Ml - m\ - 3 V )\~\U, m 2 , M 2 ) 
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+ S(m, 2 ,0,m|)| --ml{Ml-m]-s„)\- 1 {t„,mlMl)^ 

+ 5(0,TO 2 ,M|){g 

+ ^ (** - 4M 2 )(M 2 - to 2 + tOA- 1 ^, to 2 , M 2 ) } 

+ B(m 2 ,0,M 2 K ) ^2m$M 2 (s n -4M 2 )\- 2 (t n ,m 2 ,M 2 ) 

+ ^ % 2 (^ - 4M 2 )(M| - to 2 - s w ) 2 \~ 2 (t n , to 2 , M 2 ) | 

+ B(t n , to 2 , M 2 ){ + i [^(M 2 - 2m 2 - s w ) 

- {Ml-m 2 ){M 2 K -s«)]\-\t^mlMl)} 
+ B(t.,M 2 ,M|){ -g 

- ^ (** - 4M 2 )(M 2 - to 2 + ^A- 1 ^, to 2 , M 2 ) 

+ i m 2 (s^ - 4M 2 ) [-^(3M 2 - to 2 + f w ) 

+ (M 2 - m 2 )(5M 2 + to 2 - U)} \~ 2 (t n , to 2 , M 2 ) } 
+ C(M 2 ,t w ,mlO,M 2 ,M 2 K )x 

^m 2 Ml{M 2 K -m 2 )\-\t^mlMl) 

+ ? m 2 M 2 (s„ - AM 2 )(M 2 K - m 2 )(M 2 K - to 2 - S .)A- 2 (t., to 2 , M 2 ) 
+ C(to 2 , 0, m 2 , 0, to 2 , M|) x 

| _ I m 2 (M 2 _ m 2 )(M 2 + m 2 _ QX- 1 ^, TO 2 , M 2 ) | 

+ C{t n ,t n ,0,m 2 ,M 2 ,M 2 K ) x 

{ " ^ - ™?)(^ - ™f + t v )\-\t n , TO 2 , M 2 ) 
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2 

+ -^(s n - 4M 2 )(5M 2 + m 2 - t n )(M 2 n - m 2 - U)\-\t n , m 2 , M 2 ) 

2 

+ (s. - ±Ml)(M 2 K - m 2 - S7T )(3M 2 - 3m 2 - QX' 1 ^ m 2 , M 2 ) 



+ g(M 2 -m 2 + o}, 



= - — (24^ + 24fT 2 + 8K 5 + 8K 6 - 36K 12 + 12X 1 + 9X 6 ) 

lo 

I I 1 ( TTL 2 Ul 2 



+ B{MlQ,Ml)l-m 2 Ml\-\U,mlMl) 
- (M 2 + m 2 - t n )(M 2 - m 2 + QX' 1 ^, m 2 , M 2 ) 
| m 2 M 2 ( S?r - 4M 2 )(M| - m 2 - S7r )A- 2 (^, m 2 , M 2 ) 



+ \{M 2 K - m 2 - Stt )(M 2 k + 2m 2 - s ff )A" 1 (* w , m 2 , M 2 ) j 
+ 5(0,m 2 ,M|){i-g 

- ^ (s* ~ *M 2 )(M 2 - m 2 + t n )\~\t n , ml M 2 ) 

1 



+ B(m 2 ,0,M 2 K )^ ---2m^M 2 (^-4M 2 )A" 2 (t 7r ,m 2 ,M 2 ) 

- -L m 2 (s w - AMl){M 2 K - m 2 - s n ) 2 \~ 2 (t n , m 2 , M 2 ) 
+ B(U,m 2 ,Ml) | X -m\{Ml - m 2 )A^, m 2 , M 2 ) 

- 1 ^(2M 2 - 3m 2 - 2s w )A~ 1 (^, m 2 , M 2 ) J 
+ B(U,MlM 2 K ) 



mj 



(36) 



4** 
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2 
777 

+ ^ (a* - 4M 2 )(M 2 - m 2 + ^A" 1 ^, m, 2 , M 2 ) 
+ i m 2 (s, - 4M 2 ) [U3M 2 - m 2 + t n ) 

- (M 2 - m 2 )(5M 2 + m 2 - ^)] A" 2 ^, m 2 , M 2 ) J 
+ C(M 2 ,^,m 2 ,0,M 2 ,M|) x 

| -m^M^-m^A-^m^M 2 ) 

- ^ 2 M 2 ( S7r - 4M 2 )(M| - m 2 )(M| - m 2 - s^(U, m 2 , M 2 ) 
+ C(m|,0,m|,0,m|,M|) x 

| i m 2 (M 2 - m 2 )(M 2 + m 2 - f^A" 1 ^, m 2 , M 2 ) | 
+ C(^,^,0,m|, M 2 ,M|)x 

{ -^W-^ + O 
+ 1 m^M 2 - - m 2 )(M 2 - m 2 + t n )\-\t n , m 2 , M 2 ) 

2 

- ^ (a* - 4M 2 )(5M 2 + m 2 - ^)(M 2 - m 2 - ^A" 1 ^, m 2 , M 2 ) 
+ 1 m?(s ff - 4M 2 )(M 2 - m 2 - s^A" 1 ^, m 2 , M 2 ) 



3m 



2 



( Sw - 4M 2 )(M£ - m 2 - S7r )(M 2 - m\ + i 7r )A _1 (t 7r , m|, M 2 ) 

lot,,- 

(M| - m 2 - S7r )C(m 2 , i w , M 2 , m 2 , m 2 , M 2 ) . (37) 



4 The non photonic contribution 
4.1 One-point functions 

Let P denotes a pion, n, or a kaon, if, and A P the difference, 

A P = Mp± - M 2 • (38) 
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We shall expand the one-point function, 



A{M 2 ) = -i/i 



4-D 



d D l 



(2vr)^ I 2 - M 2 p0 



2 ' 



(39) 



to leading order in Isospin breaking. 

In Dimensional Regularization, the preceding integral reads, 



A(M 2 ) = M 2 p0 -2 A 



16tt 2 



In 



M 2 p0 
V 2 



By 



, this is equivalent to, 
A(M 2 p0 ) = -2A(M|± - A P ) 

{M 2 P± - A P ) In 



16n- 



M 2 p± 
V 2 



1 



Ap 

M 2 



Expanding to first order in Ap, we obtain the splitting between strong and 
electromagnetic interactions in one-point functions, 



A(M 2 o) = A(M 2 P± ) + 

4.2 Two-point functions 

The loop integral, 



16vr 2 M 



2 ml) 



A, 



(40) 



B(p ls 7no,mi) = -in 



4-D 



d U l 



(2tt) d (l 2 -m 2 )[( Pl + l) 2 -m 2 } 



21 ' 



(41) 



is function of three scalars, p 2 , TOq, and m 2 . In order to obtain Isospin 
breaking corrections generated from (}4~T|) we shall expand B(p 2 + 5, wig + 
S , m\ + 5i) to first order in 5, <5 , and Si, where these quantities are leading 
order in Isospin breaking, 



8, 8 , 8i, = 0(a, m d -m u ). 
In Dimensional Regularization, 



(42) 



16vr 2 



4-L> 



ln(47r/i 2 )+r'(l) 
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1 f 1 

— - / dx In [xml + (1 — x)m\ — x{\ — x)pf\ . 
Jo 



One then has, 

B(pl + 5, m 2 + 5 ,m? + 5i) 



16tt 2 



4-D 



+ ln(4vr/i 2 ) + T'(l) 



1 Z" 1 

167^ y dx in + ^ + ^ ~ + ^ ~ ^ ~ + ^ ' 



Expanding to first order in 5, So, and S±, the preceding equation takes the 
form, 

B(pl + 5,ml + 5 ,ml + 5i) = B(p\,ml,m 2 1 ) 
dx 



16tt 2 
1 

16^2 



dx 



dx 



X 


OTq 


+ (1 


— x)m\ — 
1 — X 


x(l 


— x)p\ 


XTTIq 


+ (1 


— x)m\ — 
—x(l — X 


x(l 

) 


— x)p\ 


OTq 


+ (1 


— x)m\ — 




— x)p\ 



■ s 1 

■ 5. 



If we denote by r the generic integral, 

T(pl,ml,ml) = / dx j 

Jo xm o 



+ (1 — x)m\ — x{\ — x)p\ 



2 ' 



(43) 



then, the splitting between strong and electromagnetic interactions in two- 
point functions is easily obtained from the following compact formula, 



B{p\ + 5, ml + So, ml + 5,) = B{p{, m 2 , m\ 



32tt 2 pI [ \mf 



ln ( — I ) + (Pi + m i _ m o) T (Pi, m l, m l) 



mi 



ln ( — 2 _ (Pi ~~ m l + m o) r (Pl 5 m 0> m l) 



32tt 2 p? 

^{ 2 p? + K-^)lnQ) 



5o 
Si 



32tt 2 p 

+ [(m 2 - ml) 2 - p\(m\ + m 2 )] t{j>\, m 2 ,, m 2 ) }> 5 . 
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(44) 



As an application, consider the two-point function, B(pi — p, M w o, M K o), 
selected from the t-channel contribution to A°~. The following replacements 
in (JIIJ), 

5 — > -— {Ml + m 2 - t 7r )A 7r - 1 A A - 



8n 



-A, 



2 



lead to the expression, 



B( Pl -p,M w0 ,M K o) = B(t n ,Ml ± ,M 2 K± ] 



+ 
+ 
+ 



A. 



327T ijj- I 5 jy 



1 + — (M 2 + m 2 - t.)(M| - Ml) 



111 



M| - M 2 + U ~ - {M\ + Ml ) 
(M 2 + m 2 - ^)(M| - Ml) 2 



T(t n ,Ml,M 2 K ] 
Ml 



+ ^k\ 1+ ^- M "- 2t ^\M 2 K 



- (3M 2 K -Ml- 2Q 



r(U,Ml,M 2 K/ 



(45) 



4.3 Isospin limit 



ui = o 

U! = 



24F 2 1 16vr 2 



4(2M 2 K± + 4M 2 ± - s w ) 
3 



+(2M^± + 4M 2 ± -U)-- (M 2 ± - M 2 K± )(Ml ± + M 2 K± ) 
+ (QM 2 K± + £ (M 2 - M| ± )(M 2 + M| ± ; 



+ 48 



2(M 2 ± + M 2 K± - U)L 3 + 2(M 2 ± + 2M 2 K± )L A - m 2 e Lc 



(46) 
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+ a{m;±) 
+ ml) 

+ A{M 2 K± ) 



5 - £ Ml ± - UMl ± - M 2 K± ) 
_ 3 _ 6 (5M 2 _ 6M ^ ±) + 18 (M 2 _ M ^ ±)5 

-2 + ^(M 2 ± -A4 ± ) 



+ r K - M k±) + | " Ml ± f - ^ {Ml - M 2 K± f 

+ 4(4M 2 ± - S7r )B( S7r ,M 2 ± ,M 2 ± ) 
+ 2{AM 2 K± -s- K )B{s 7n M 2 K± ,M 2 K± ) 



+ B(U,M 2 ± ,M 2 K± ) 



-G(M 2 ± -M 2 K± +U) 



+r( M ^- M x ± ) 2 + ^« ± -M^ ± ) 3 



+ B{t^M 2 q ,M 2 K± ) 



6(M 2 - 3M|± + f w ) 



+1 ( M J - M| ± )(5M 2 - 3M| ± ) - i? ( M 2 - Mj^ 



The e-terms 



+ 



+ 



-3 
1 1 



24F 2 16tt 2 
15 



2(6M 2 + 28M 2 + 20M 2 - 9U) 
9 



+ - (M 2 - M 2 K )(M 2 + M 2 K ) - - (M 2 - M 2 K )(M 2 + M 2 K ) 
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+ "w7 
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r L 



2(M| + 5M 2 -2^-^)L 3 



+6(M 2 + 2M 2 )L 4 + 48(M£ - M 2 )(3L 7 + L 8 ) - 3m 2 L 9 



8F 2 



A(M 2 ) 



32M 2 



M 2 - M 2 



2 10 
-15 + -M 2 --(M 2 -M|.) 2 
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+ 3A(M%) 
+ 2A{M 2 K ) 
+3 + 



U? = 



3 + A (3M 2 - AMI) - I {Ml - M 2 K ) 2 
16M 2 



M 2 - M 2 

TV T) 



i (5M 2 K - AM 2 ) - A (M 2 - M 2 K ) 2 



-A(2M 2 -3M|) + ^(M 2 -M 2 

12(2M 2 - S7r )£( S7r ,M 2 ,M 2 ) 
+ 2(4M|-3 S7r )S( S7r ,M|,M|) 
+ 4(3M 2 + M 2 - 3s v )B(s n , M^M 2 ) 



2 \2 



+ 2B(t w ,M 2 ,M 2 K 
1 



7M| + 13M 2 - 9t w 



2 \3 



+_ (M 2 _ M £)(5M| + M 2 ) - - (M 2 - A^) 



+ 2B(t w ,MlM 2 K ) 3M 2 - 5M 2 K + 3t w 

A{Ml) 
A(M 2 V ) 



AF 2 



1 + -M 2 K 

b-rr 



r (3M 2 



AM 



K / 



^ M| + 1 (3M 2 - 4M|) 



A(M|) 



^M 2 K {Ml-M 2 K )B{t^MlM 2 K ) 
+ 2B(t n ,MlM 2 K )x 

3(M 2 - M 2 K ) - i (M 2 - M|)(3M 2 - 4M|) 
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4.5 The A .--terms 
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1 
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1 


256tt 2 F 2 





(4M| + 20M 2 - 3s u ) 



4(M 2 + m 2 -0 



+ r (M 2 + M 2 )A«-1(M 2 + 6M| + M 2 )^ 



+ - [(5M 2 - 3M 2 )A^ - 5(3M 2 - M|)A^] A tK 

7T 
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+| A^ + i (M| - 2Ml) A lK - ^ A lK A\ K 



tl 



_6 



1 1 



8F 2 Sw 



1 1 A 



+^ (6M 2 - 7M 2 K )A eK + - AffArfA^ 



+ 



f2 

7T 

1 1 



24F 2 s n 
9 



9 3 



6 



+- (M 2 - 2M£ + M 2 )A, X + - (A nK + A vK )A m A nK 



iln 



- M|s, + - M 2 K A nKSlT 



2567r 2 F 2 s w 

-1 (M 2 - + 3M 2 ) A m A nK + | A W A^ + 1 A w A 4 rf 



+ 



5127r 2 F 2 s w 



lln 



- (3M 2 - M|)A«A^ 



-1 (3M 2 - 2M 2 K )A m Al K - | A^A^A^ 
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256tt 2 F 2 s n 



2M 2 K s n - - (M 2 K + 3M 2 )Z nK A eK 
+i (3M 2 K + 5M 2 )A eK Al K - I M 2 K A\ K s % 
+i (3M 2 K - M 2 )A eK Al K - I A eK Al K 



1 



1 



— r(f w , M^, M£)A^(A^ - E^) 



x 



512tt 2 F 2 s w 

| (3M^ 2 - M|) - 1 (3M* - 2M|) - A A^A 
J_ 1 (5M 2 - 2^)5^, M 2 , M 2 ) 

O-Tq S,r 



-L±(M 2 K - S77 )B(s„,M 2 K ,M 2 K ) 



1 1 



B(U,M 2 ,M 2 K ) 



1 

^7T 



(5Ml - 2M 2 )s„ 



SF 2 s n 

r 7T r 7T 



1 1 



16F 2 s n 



-B(t w ,M^M 2 K ) 



^M 2 S7T -^A 2 vK s n 



+ ^ (3M 2 - 2M 2 K )A eK A vK + 1 A^A^A^ 



f 2 

8 



— 12(M 2 — \.M\ + M 2 — 3s,r) 

6(2M 2 K - 3M 2 - m 2 )A wK + (M 2 K + 5M 2 )s n 
-1%M 2 K {3M 2 - %M 2 K - m 2 ) + 3{M 2 + M 2 K )s n 
-6M 2 {7M 2 + 2M 2 K + m 2 ) 



768vr 2 F r 
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+ 



+ 



+ 



tl 



4A 



+3M 2 K (7M 2 K -M 2 + 3m 2 )A wK 
-3M 2 (2M 2 K - 2M 2 + m 2 )A wK 
-9Ml(10M 2 - 5M 2 - m 2 )A vK 

+3M 2 (66M 2 -81M 2 K - 13m 2 ) A vK 
-^\(M 2 K + 2M 2 + m 2 )Al K 
+3(6M 2 -9M|-m 2 )A^ 
-^-A(M 2 ) 

24F 2 s w i jvt; i v 

~ (2M 2 K - 3M 2 )M 2 K + 1 (3M| - 2M 2 + m 2 )M 2 

- i- (5M 2 - 6M| - s ff ) - 72 A ^ s - 
-| (ISM 2 - \1M 2 K - 5m 2 )M 2 



-1 (20M 2 - 21M| - 3m 2 )M| 
12 

--(M| + 2M 2 + m 2 )A^ 



1 



-A(M 2 K )x 



48F 2 s n 

- ^ (6M 2 - 12M 2 - 7 S?r ) + £ (4M 2 - 6M 2 K - s w ) 
- {5A nK s w - 6Z eK M 2 - 15A wK M 2 K + 3A £n M 2 K ) 

7T 
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+ - A vK s n - - (3M 2 K + 2m 2 )M 2 
18 

+ -(5M 2 K + m 2 -3M 2 )M 2 K 
-^(M 2 K + 2M 2 + m 2 )Al K 
-^(QM 2 -9M 2 K -m 2 )A 2 )K 

+ 76SW^ h '(f){- 2(3M -- 3M - +2S - ) 
+i [2(2*4 + 3M*)s„ - 3(3*4 + m])A, K ] 

+i [2M 2 K s n + 3(M 2 K - 4M 2 - ^ 2 )A^] A nK 

+ |-(M| + 2M 2 + m 2 )A 4 ^} 
1 / M 2 \ 



-A 3 



~ 192^2^ i^ M l ~ s -jr)r(s n , M 2 , Ml) 

l 2(33W - 2 S?r ) + 1 (Ml + 2M 2 + m 2 )A^ 

~ [2Kks* - 3M 2 A eK + 3M 2 K (5M 2 K + 3m 2 )} A 2 nK 

-1 [3M 2 (A« + 4M 2 ) + M 2 K {21M 2 K + 9m 2 - 2 S?r )] A nK 

7T 

--^ [3M 2 (A« - 2 S?r ) + Ml(21M 2 + 3m 2 - As w )] ) 
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+ 



512tt 2 F 2 



r(^,M 2 ,M|)(A 2 ;i ,-£^) x 



r + 1 K - 3M|) - |- (5M 2 - 3M|)A^ + 



+ 
+ 



3F 2 



4 

B(s n ,M 2 ,M 2 ) 



1 1 



-B(t n ,M 2 ,M 2 K ) 



2s, 



+ 



24F 2 S?r 
+ r [6A^-(7M|-4M^) S J 

(15M| - 6M 2 + 3m 2 - 7s n ) A 2 wK 



--(M 2 K + 2M 2 + m 2 )Al K 



T [2(5M 2 - 3M|)A^ - (2M 2 - M 2 K )s,] 
- M 2 (18M 2 - 9M 2 - 5m 2 ) 



-A^-s* + 3M|-(10M 2 - 5M 2 K - m 2 ) 



12 



+ -(M| + 2M 2 + m 2 )A^ 



In the preceding expressions, 



A £P 


= "i/ 


-M 2 P , 


S^p = 


Apq 


= M 2 P 


~M 2 Q , 


SpQ = 



M 2 + M 2 . 



4.6 The A^-terms 



768tt 2 F 2 



- 16 



M 2 + 2M 2 K 
M 2 - M 2 
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+ 12 + - (3M 2 - 6M 2 K - Ml) 



+ ^ [2(M| - 3M 2 ) A„ K - 9AJ„] ~^l K ) 



24F 2 tl y nJ 



7T 
1 1 



24F 2 M 2 
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384n 2 F 2 
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3 + -(4M 2 K -Ml) 



+1 M 2 k Kk -h^K- 3Ml)Al K + I A^ 
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512tt 2 F 2 \M 2 k 



M 2 
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- (3M 2 - M|) 



1(3M 2 -2M|)A^ + J-Aj 



tl 
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384vr 2 F 2 



t(U,MI,M 2 k ) 



y-3 ^ 

- (7M 2 - Ml - 3t w ) 
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A 4 - - 
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24F 2 
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IT 
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4+ (4M 2 -5M^)--A^ 
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768tt 2 F 2 



6 + - (34M| + 9M 2 - 3M 2 ) 



+- [(17M 2 K - 3M 2 )A„ K - 3(7M 2 + 3M 2 K )A r)K ] 

TV 
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8F 2 



16F 2 
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ml) 
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4M 2 

7T 

Ml 
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M| 
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384tt 2 F 2 



768tt 2 F 2 



tl 
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7T 
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7T 
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+ 



+ 



7T l TT 

+ 512^F 2 T(t ^ M '' (~ 2M ^ + ^ + ^ A ^ 
2 + A (M 2 - 3iW|) - | (5M 2 - 3M|)A^ + 



x 



^A 3 " 

tl ^ K 



6F 2 



8F 2 



B{s^M 2 K ,M 2 K ) 
B(t n ,M 2 ,M 2 K ) x 



3 2 S#2 A 5 a2 2 a3 



-2 + — (2M| - ill 



16F 2 



B(U,M 2 ,M 2 K )x 



2 - f (2M 2 - 5Ml) - 1 E^A^c + ^ A 3 ^ 
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(55) 



5 Results 



5.1 Input 

The numerical values of the physical parameters must be fixed through ex- 
perimental input. However, this input may not necessarily consist of direct 
measurements of the renormalized parameters; it may be obtained from any 
suitable set of experimental results. In practice one uses those experiments 
which have the highest experimental accuracy and theoretical reliability. This 
criterion is certainly fulfilled for the following set of parameters whose nu- 
merical values are taken from [TU] : 

• the fine structure constant, 

a = 1/137.03599976(50) , 
corresponding to the classical electron charge e = y/4na, 

• the masses of the charged leptons, 

m e = 0.510998902(21) MeV , = 105.658357(5) MeV , 
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• the Fermi constant 



Gp = 1.16639(1) • 1(T 5 GeV 



-2 



which is directly related to the muon lifetime, 



• the Cabibbo-Kobayashi-Maskawa quark-mixing matrix element, 



its 



0.2196 ±0.0026, 



coming from the analysis of K e3 decays, 



• the masses of the light mesons, 



M,± = 139.57018(35) MeV , M w o = 134.9766(6) MeV , 

M K ± = 493.677 ± 0.016 MeV , M K o = 497.672 ± 0.031 MeV , 
M v = 547.30 ± 0.12 MeV , M p = 771.1 ± 0.9 MeV , 



• the charged light mesons decay constants, 

F n ± = 92.419 ±0.325 MeV, F K ± = 112.996 ± 1.301 MeV , 
coming from the analysis of ir^ and K^ 2 decays, respectively. 

Let us consider the parameters, M n , Mk and e, related to light quark 
masses. Since and Mk figure in our expressions only at next-to-leading 
order, it is completely safe to replace them by their leading order expressions. 
In fact, the quantity M„ will be identified to the neutral pion mass, 



7T 



134.9766(6) MeV 



while M K will be replaced by, 




- (M K± + M K0 + Ml, - Ml ± ) 




495.042 ± 0.034 MeV . 



e = 



(1.061 ±0.083) ■ 10 



-2 
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extracted from the mass splitting in the baryon octet. 

We will turn now to the determination of low-energy constants in the 
strong sector. Following [T2], these constants will be evaluated at one- loop 
accuracy, that is, by fitting experimental measurements of the concerned 
observables to their ChPT expressions at next-to-leading order. Note that 
all of our expressions will be evaluated at the scale \i equal to the rho mass. 
The form factors are sensitive to variations of the low-energy constants, 
L\, L2 and L3. By fitting experimental results on form factors |T3] to 
their ChPT expressions at next-to-leading chiral order we obtain [14j, 

L\ = (0.46 ± 0.24) • 10~ 3 , 

V 2 = (1.49 ± 0.23) • 10~ 3 , L r 3 = (-3.18 ± 0.85) ■ 10~ 3 . 

The constant L5 can be extracted from the ratio of the kaon to the pion 
decay constant in the isospin limit [T2] . 

Fk± - 1+ 4 fM 2 MM 77 -I 5M ^ ln M - 2± 
~ 1 + Ff ± {Mk± ~ M * ±)L * + 128tt 2 F 2 ± ln ~f 

M 2 K l M 2 K± 3M 2 m\ 

64tt 2 F 2 ± 11 fi 2 128vr 2 F 2 ± n /z 2 ' 

and reads, 

L\ = (1.49 ±0.14) • 10~ 3 . 

Having L 5 , it is easy to determine L 8 from the quantity accounting for 
5*7(3) breaking [TJ], 

Am = -^{M 2 K± -Ml ± ){2Ll-LD 

M 2 . , M 2 M 2 , M 2 
- ln — — A - — ln — - 

32vr 2 F 2 ± /i 2 32vr 2 F 2 ± /1 2 ' 

and which value reads [TTj, 

A M = 0.065 ±0.065. 

The result is, 

Lg = (1.02 ±0.22) • 10~ 3 . 



28 



The constant L 7 is obtained from L 5 and L s with the help of the isospin limit 
quantity, 

A GMO = (4M£± - Ml ± - 3M*)/(M* - M* ± ) = 0.2027(15) , (56) 
by matching its value to the ChPT expression at next-to- leading order [T2*] . 

Agmo = -^r(M*-M* ± )(l2L 7 + 6Ll-L r 5 ) 

2 / M 4 ± M 2 ± 
+ — — In — 

M 2 -M* ± \32n 2 F* ± /i 2 

8tt 2 F 2 ± /i 2 32tt 2 F 2 ± /x 2 J 

We obtain for L7 the value, 

L 7 = (-0.44 ±0.12) • 10~ 3 . 

The constant L 9 is fixed from the electromagnetic charge radius of the pion [T5] . 

Lg = (5.5 ±0.2) • 10" 3 . 

Finally, it is difficult to fix the constants L 4 and L 6 by direct experimental 
determination. These constants are suppressed by the Okubo-Zweig-Iizuka 
(OZI) rule and measure the amount by which m s affects the values of the 
order parameters F and (qq). The constant L 4 was derived from Roy and 
Steiner equations for S- and P- waves of it — K scattering amplitude JHj , 

LI = (0.53 ±0.39) • 10~ 3 . 

The constant L$ has been obtained from a chiral sum rule |17j . 

LI = (0.4 ±0.2) • 10" 3 . 

To close the discussion about the strong sector we have to fix the param- 
eter F . At leading chiral order this parameter is given by the pseudoscalar 
decay constants, F n , F^ or F v . One can then see the latter as the "renormal- 
ized" quantities corresponding to the "bare" quantity F and thus replace 
it by one of them after accounting for next-to-leading order contributions. 
But the main question is which expression for the decay constants to use 
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(57) 



especially that the difference between their numerical values is relatively big. 
For instance, the expressions for the pion and kaon decay constants at next- 
to-leading order are given in the isospin limit by 



F K± 
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[ 1+ F 

1- 7 


M 2 ± 

TV 1 


16tt 2 F 2 ± 




L 7 


3M 2 ± 

TV 
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iir 2 F 2 



{Ml ± +2M 2 K± )U A + 



AM 2 , 



F 2 



jr 



In 



M 2 



M 2 K± 
32n 2 F 2 ± 



In 



M 



K± 



(M 2 ± +2M 2 K± )Ll + 



4M|± 
F 2 



In 



M 2 



3M 2 K± 
6An 2 F 2 , 



n 2 



3M 2 



128vr 2 F 2 ± 

7T -J- 



M 2 
In— ^ 

/i 2 



and L5, we 



Taking as input the aforecited values for M^i , Mk± , , , L\ 
obtain for Fq the central values, Fq = 67.53 MeV and Fq = 57.40 MeV from 
F^± and F K ±, respectively. If, for comparison, we take for L 4 its large- iV c 
estimate, the central values modify to Fq = 79.16 MeV and Fq = 71.62 MeV 
from F n ± and F K ±, respectively. This amounts for a 15% to 20% deviation 
for the value of F . In our calculation we will use for F the two values given 
by the bounds of the following inequality, 

57.40 < F < 67.53, 

and give the difference between the two obtained results as an error on the 
final result. 

In the electroweak sector it is quasi impossible to have an experimental 
determination of the low-energy constants due to the relatively big number 
of constants from one side and to the relatively small magnitude of the elec- 
troweak effects from the other side. We will use for the constants Ki in the 
mesonic sector the following central values obtained by means of resonance 
saturation jTH], 



K[ = -6.4 



-6.4 



io~ 3 
io~ 3 



K\ 
Kl 



-3.1 ■ IO- 3 
19.9 ■ 10- 3 , 



Kl 



Kl 







K 



10 



0. 



K 
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: 6.4- 10~ 3 , 

: 8.6- 10~ 3 , 
-9.2 ■ 10~ 3 



with an error of ±6.3 • 10 -3 assigned to each of them coming from naive 
dimensional analysis. The latter will also be used to fix the bounds on low- 
energy constants in the electroweak leptonic sector, 

IXI < 6.3 • 10~ 3 , 
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since these constants have not been yet determined. 



5.2 The / form factor 

In what follows we will refer to the inequality, 

4M 2 ± <s n < (M K ± - m t f , (58) 

as the allowed kinematical region. The first term in the partial wave ex- 
pansion for / form factor is infrared finite. It contains however singular 
(Coulomb) terms for, 

s n = [M K - mtf + 2m t {M K t 2M n - m t ) . (59) 

As can be easily seen, the singularity is outside the allowed kinematical region 
for me ^ and approaches the upper bound from the right when me tends to 
zero. Therefore, there is no apparent reason for subtracting Coulomb terms 
in the case of non- vanishing lepton mass. In order to see the impact of such 
terms on the whole correction, let us consider the following imaginary part, 

3 f f 1/2 

lmU f ( S7r ) = — 4= K - 2M " 



32ttF 2 v/3 v * nJ \ s 



48ttF 2 



+^ ^ (t, + Ml - mJ)\-^(U, ml Ml) . (60) 

OZ7T t w 

The plot of the preceding expression as a function of is given by figure El 
It is easy to see that e 2 (singular) -terms are almost negligible with respect 
to A T or e-terms. 



5.3 The g form factor 

Unlike the / form factor, the g form factor is infrared divergent. We have 
shown in pQ that this divergence is cancelled at the level of differential decay 
rate by the one coming from real soft photon emission. In experiments, 
one has to measure modules and phases for form factors. Therefore, a sub- 
traction of the infrared divergence should be applied at the level of form 
factors. The trouble is that the subtraction is not unique. A possible choice 
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corresponds to a minimal subtraction and consists on dropping out the In m 7 
term. Another possible choice which we qualify by a reasonable subtraction 
consists on treating / and g form factors on a equal footing. While the / 
form factor is infrared finite, the infrared divergence in g form factor comes 
from wave function renormalization of external charged particles and from 
virtual photon exchange. The latter contribution is generated from the Cq 
function, 

Co(-pi,P2,m 1 ,mi,M 7r ) , (61) 

expressed by formula (195) in the appendix of reference pQ. In the reasonable 
subtraction scheme, one drops out the In m 7 term coming from wave function 
renormalization and the full contribution of the Co function. Formally, one 
introduces a subtraction parameter, £, which equals 1 in the minimal sub- 
traction scheme and vanishes in the reasonable one. Having this, we define 
the subtracted real part, 

g P {s v ,Z) = l + ReU 9 (s n ) 

e 2 U-Ml-m 2 

:f x 



2 



8tt 2 ^/ tn - (m e + M n ) 2 ^t n - (m e - M n f 

y/t n - (m e - M^f + y/t n - (m e + M n ) 2 

In — — In m 

y/U - (m e - M„) 2 - y/U - (m e + ikQ 2 

+2e 2 (t n - Ml-m 2 e) x 

(1 - Re C(ml t n , Ml m 2 , m 2 , M 2 ) . (62) 



Finally, from the imaginary part, 
ImC/»( S .) = «}(«„) 



32nF 2 V s» 



1/2 



1 X-y^mlM^x 



327rt w 

[m 2 e (5U + M 2 - m 2 ) + At w (M 2 - Q] 

-2e 2 (t n - Ml - m 2 ) Im C(m|, t n , M 2 , m 2 m 2 , M 2 ) , (63) 
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where, 



o 

we define the subtracted phase as, 
5 P (s„,Z) = lmU 9 ( SlT ) 



+ ^ % - M l - ^ 2 )A- 1/2 (t., m £ 2 , M*) i him 2 
+2e 2 (^-M^-m^) x 

(1 - Im C(ml t n , Ml ml, m 2 e , M 2 n ) . (65) 
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6 Conclusion 

In this work we proposed a possible splitting between strong and electromag- 
netic interactions in decay form factors. The technique was applied to 
the decay of neutral kaon, K° — ► 7r°7r~£ + z/^. It consists on working at the 
production threshold for the lepton pair, S£ = mj. The latter assumption 
simplifies significantly the splitting by allowing a partial wave expansion of 
form factors with exactly the same structure as in the pure strong theory. 
This constitutes a good approximation as long as the dependence of form 
factors on se remains linear; the slope poor. 

The interest in the present process is at first theoretical. In fact, the 
partial wave expansion of form factors involves the P-wave iso-vector 
phase shift, 5l(s n ), which can be related to txtx scattering lengths via Roy 
equations. In their turn, scattering lengths are sensitive to the way Chiral 
symmetry is spontaneously broken. Consequently, a theoretical study of the 
process in question including all possible contributions is imperative. We 
gave here the first analytic and numerical evaluation of the Isospin breaking 
contribution. This would allow the extraction of Sl(s n ) from the experimental 
measurement of form factors. 

We started with the evaluation of the first term in the partial wave expan- 
sion for / form factor. This term vanishes in the absence of Isospin breaking 
and is free from Infrared divergences in its presence. Motivated by these two 
features, we studied the sensitivity of Isospin breaking correction to varia- 
tions of F and rri£. This was achieved by plotting the graph of the correction 
as a function of s w for two values of F in figured] and for = m e) m^ in 
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figure |21 We then compared in figure El the relative size for the different 
contributions to the correction coming from virtual photons, 0(e 2 ), mass 
square difference between charged and neutral mesons, O(Z e 2 ), and mass 
difference between up and down quarks, 0{md — m u ). 

We pursued with the evaluation of the first term in the partial wave 
expansion for g form factor. The comparison between the size of Isospin 
breaking correction to the real part of the term in question and the one-loop 
level correction to the same quantity and in the absence of Isospin breaking 
was made in figure |U Oira^ — m u ) and 0(a) contributions to the preceding 
correction were compared in figure Finally, Isospin breaking correction to 
the P-wave iso- vector titc phase shift was plotted in figure |H1 

Our results are of great utility for the interpretation of the outgoing data 
from the KTeV experiment at FNAL. 
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For the following integral, we shall distinguish between two 
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The lepton is a muon. 
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A. 2 r-integrals 

These integrals appeared while splitting strong and electromagnetic parts 
in two-point functions. Their definition is given by equation (|43|). We are 
interested in the following particular r-integrals. 
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A. 3 C- integrals 



These are scalar three-point functions whose definition and expressions were 
given in the appendix of reference [T|. In what follows, we sketch some of 
the particular cases which we need for the numerical evaluation of Isospin 
breaking corrections. 
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Re U f (s„) 




Figure 1: The real part of the first term in the partial wave expansion for 
/ form factor under the assumption, sg = mj = w? e . The error band comes 
exclusively from the uncertainty in the determination of low-energy constants 
and has been developed in quadrature. 
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Figure 2: The real part of the first term in the partial wave expansion 
for / form factor under the assumptions, se = mj, F = 67.53 MeV. The 
error band comes exclusively from the uncertainty in the determination of 
low-energy constants and has been developed in quadrature. 
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Figure 3: The imaginary part (in Radian) of the first term in the partial 
wave expansion for / form factor under the assumptions, se = mj = m%, 
F = F w = 92.419 MeV. 
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Re U 9 (s n ) 




Figure 4: Radiative correction to the real part of the first term in the partial 
wave expansion for g form factor under the assumptions, Se = mj = ml, 
F = 67.53 MeV. The infrared divergence has been removed applying a 
minimal, £ = 1, as well as a reasonable, £ = 0, subtraction scheme. Error 
bands come exclusively from the uncertainty in the determination of low- 
energy constants and have been developed in quadrature. 
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Figure 5: Isospin breaking correction to the real part of the first term in the 
partial wave expansion for g form factor under the assumptions, se = mj = 
?T7.g, Fq = 67.53 MeV. The infrared divergence has been removed applying a 
reasonable, £ = 0, subtraction scheme. Error bands come exclusively from 
the uncertainty in the determination of low-energy constants and have been 
developed in quadrature. 
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Figure 6: The imaginary part (in Radian) of the first term in the partial 
wave expansion for g form factor under the assumptions, Si = mj = m^, 
F = F w = 92.419 MeV. The infrared divergence has been removed applying 
a minimal, £ = 1, as well as a reasonable, £ = 0, subtraction scheme. 
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